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The local viscous-inviscid interaction field generated by a wall temperature jump 
on a flat plate in supersonic flow and on the windside of a Reusable Launch Vehicle in 
hypersonic flow is studied in detail by both a Navier-Stokes numerical code and an an- 
alytical triple-deck model. Treatment of the rapid heat transfer changes both upstream 
and downstream of the jump is included. Closed form relationships derived from the 
triple-deck theory are presented. The analytically predicted pressure and heating varia- 
tions including upstream influence are found to be in generally good agreement with the 
Computational Fluid Dynamic (CFD) predictions. These analyses not only clarify the 
interactive physics involved but also are useful in preliminary design of thermal protec- 
tion systems and as an insertable module to improve CFD code efficiency when applied 
to such small-scale interaction problems. The analyses only require conditions at the 
wall and boundary-layer edge which are easily extracted from a baseline, constant wall 
temperature, CFD solution. 


Nomenclature 

C Chapman-Rubesin parameter, pT e /p e T 

C p specific heat 

Ci _6 constants in analytical solutions 

H total enthalpy 

L reference length to discontinuity 

M - Mach number 

n distance along normal to wall, m 

p pressure, N/m 2 

Pr Prandtl number 

q w wall heat transfer rate, W /m 2 

ReL Reynolds number, p e V e L/p e 

T temperature, K 

V velocity, m/s 

x,y,z Cartesian coordinates, m 

x,y coordinates, origin at discontinuity 

x j streamwise coordinate of jump 

(3 (M 2 - l) 1 / 2 

<5 boundary-layer thickness 

e scaling parameter, 

AT W wall temperature jump, T1+ — T~, K 

k upstream influence parameter = 0.827 

A Blasius solution constant = 0.332 

p coefficient of viscosity 

ui viscosity temperature-dependence exponent, 

(p cc T u ) 
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p density, kg/m 3 

Subscripts 

AD adiabatic wall conditions 

e conditions at boundary layer edge 

0 undisturbed boundary layer ahead of 

interaction zone 

REF based on Reference Temperature 

3^ = .50 + .039M 2 + .5^ 
w wall surface conditions 

oo freestream conditions 

+,- conditions downstream and upstream 

of jump, respectively 

Superscripts 

() nondimensionalized by triple-deck theory 

Introduction 

Surface property discontinuities on a hyper- 
sonic vehicle are most often associated with the junc- 
tion of two different thermal protection materials. 
They may also be associated with abrupt loss or de- 
position of surface coatings or changes in thermal con- 
ductivity and heat capacity of sensors on a test article. 
The abrupt jump in surface properties (catalytic effi- 
ciency, emissivity) in the presence of dissociated flow 
in hypersonic flight will lead to abrupt changes in heat 
transfer rate. The interactive thermophysics induced 
by such surface property jumps are remarkably similar 
to the phenomena induced by an artificially imposed 
temperature jump; consequently, insight into the cat- 
alytic jump problem may be gained through study of 
the temperature jump problem. 1 
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While' the induced heating rate associated with dis- 
continuous surface properties can be very large, the 
length scale of the interaction is usually poorly re- 
solved by conventional computational fluid dynamic 
(CFD) grids. A purely numerical analysis of the prob- 
lem can be addressed by a local grid refinement in 
the streamwise direction until the interaction is fully 
resolved. However, a better approach - one that pro- 
vides physical insight and reduces the computer time 
required for design - is to couple a local, analytic solu- 
tion of the interaction region around the discontinuity 
with a CFD solution obtained on a conventional grid. 

This paper reviews the progress to date on the com- 
putational and analytical methods used to analyze the 
catalytic and temperature jump problems as presented 
in previous papers. 1,2 The most recent analytical work 
on the temperature jump problem has revealed the 
necessity of going to a triple-deck analysis in order 
to pick up the important upstream interactions that 
were missed in a simpler integral-boundary-layer ap- 
proach. 1 Finally, it presents for the first time a test of 
the analytic, triple-deck analysis of problems of more 
general relevance to hypersonic, lifting bodies. The 
original triple-deck and computational analysis consid- 
ered only a flat plate, for which analytic solutions in 
the transformed coordinates can be obtained. In the 
more general application, a computationally derived 
coordinate - one that is matched to the baseline CFD 
solution - will be used within the triple-deck analysis to 
match the physics in the interaction region. This ap- 
proach is guided by a similar coupling procedure that 
was used with great success in problems with global 
changes in surface catalysis. 3 

The long range goal of this work is to provide an- 
alytic tools coupled with computational data sets to 
be used in the design of hypersonic vehicles in an 
Intelligent Synthesis Environment (ISE). Given a sin- 
gle, baseline CFD solution, the analytic tools allow 
a designer to investigate effects of continuous global 
changes or discontinuous local changes to thermal pro- 
tection system surface properties on heating. 

Problem Definition 

As schematically illustrated in Fig. 1 we consider an 
incoming non-adiabatic laminar boundary layer with 
supersonic external flow that encounters an arbitrary 
wall temperature jump ( T~ to T+) at some arbitrary 
streamwise station L. The incoming boundary layer 
profile is modeled by self-similar compressible flat plate 
velocity and total enthalpy profiles, while the over- 
laying disturbed external inviscid flow is assumed to 
be governed by linearized supersonic theory. Further- 
more, the Reynolds number Rei is taken to be very 
large such that the small parameter t = Re L ' is 
small compared to unity along with the boundary-layer 
thickness ratio So/L ~ e 4 . Finally, as confirmed by 
numerical results, 1 the effects of streamwise diffusion 
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Fig. 1 Schematic of interaction generated by a 
wall temperature jump. 

may be neglected even within the very short scale of 
the jump-generated interaction zone as long as the gas 
flow behaves as a continuum. 

We note that this problem is phenomenologically 
similar to the case of an abrupt jump in surface cat- 
alytic efficiency. Such a case was originally studied 4 
using integral boundary-layer theory and compared to 
computational fluid dynamic (CFD) simulation. The 
results of that study are presented here in Fig. 2. Grid 
convergence of the heating rate across the discontinu- 
ity is evident in the CFD results designated by the 
symbols. Integral boundary-layer theory, given by the 
solid line in Fig. 2, assumes that the jump in heat- 
ing and subsequent relaxation across the discontinuity 
may be predicted without consideration of upstream 
influence. Ignoring upstream influence leads to a large 
overestimate of the jump in heating rate and the re- 
laxation distance to the asymptotic heating level. An 
ad-hoc correction for upstream influence in the inte- 
gral boundary-layer theory, given by the dashed line, 
bounds the CFD results from below. 

This result suggested that a triple-deck analysis 
which implicitly includes the upstream influence would 
be a more appropriate approach for problems involving 
discontinuous jumps in boundary conditions. A triple- 
deck analysis of a supersonic flow involving a tempera- 
ture jump was chosen as a more tractable “pathfinder” 
before embarking on the more complex triple-deck 
analysis of a. nonequilibrium, catalytic jump problem. 
The derivation of this triple-deck analysis is originally 
presented in an earlier paper. 1 This paper focuses on 
the application of the theory to more general problems. 

Triple-Deck Formulation 

Details of the triple-deck formulation are available in 
an earlier paper 1 and are not repeated here. Only the 
final analytic relations describing the upstream and 
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Downstream Relations: x > 0 


p{x) = C 3 


(AT W \ r°° (t + r l /*)e- Kli 
V T w ) Jo t*/> + t*/ 3 + 1 


dt (8) 
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for large x (9) 
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ex 1 / 3 
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( 10 ) 


where 


Fig. 2 Grid refinement study of convective heat- 
ing on sphere in vicinity of discontinuous jump 
in surface catalycity and comparison to integral 
boundary-layer theory, (see Reference 4) 


C 3 = V3 (0.729)/27rP r 1/3 K 2/3 (11) 

C 4 = V3 (0.729)r(2/3)/2 7 rP r 1 / 3 « 4/3 (12) 


downstream solution for pressure and heating are pre- 
sented. The appropriately non-dimensionalized vari- 
ables with x equal to 0 at the jump as noted in Fig. 1 
axe defined: 


x = 


[x/L\\ 5 Hp 3 / 4 


e 3 C 3 J£ F (T R EF/Too) 3 /HT w /T REF r+H* 
- _ [(p - Poo)/PooV£]ff 1/2 


p = 




9 = 


Qw 

Qw,o 


(1) 

( 2 ) 
(3) 


The viscosity exponent u> is set to 0.75 herein. 
Upstream Relations: x < 0 

/at u 


p(x) = Ci y-f — i e " 


q — i - C 2 


(at u 


\ t w 


where 


Ci = (3/4) (0.729 )/Py i K 


l/3^2/3 


(4) 

(5) 

(6) 


C 2 = (3/4) (0.315) [1 + 0.196(1 - P r )] (7) 

and other referenced variables are defined in the 
Nomenclature. 


C 5 = 


V3 (0.729)r(l/3)/3 1 / 4 

2n\ 1 / i Cftg F (T w /T 00 ) 1 / 2 

Co = |V3 C 2 /2tt 


(13) 


(14) 


The integrals in Eqs. 8 and 10 are evaluated numer- 
ically. 

Coupling to CFD for General Applications 

For the flat plate test case, freestream conditions 
(subscript (,*,)) are nearly identical to boundary-layer 
edge conditions (subscript ( e )). In more general appli- 
cations we assume a locally two-dimensional flow over 
the domain in which the interaction occurs. In these 
cases edge conditions replace freestream conditions for 
all quantities used in Equations 1-14. Furthermore, 
the background reference pressure used in the defini- 
tion of p in Eq. 2 and the reference heating heating 
level q Wt0 used in Eq 3 are replaced by the locally 
varying value p e (x) and q w (x) from the baseline CFD 
solution. Use of these baseline reference values allows 
us to isolate the perturbations caused by the jump 
condition. 

Also in the general case, the reference length L at 
jump location xj required for the Reynolds number 
Re 1 is derived from the integral boundary-layer the- 
ory expression for the local, undisturbed value of q w 
in the baseline CFD solution. This approach is identi- 
cal to the algorithm used to compute effects of global 
changes in surface catalytic efficiency on heating from 
a baseline CFD solution. 3 Consequently, given a base- 
line CFD solution for the heating rate q w the reference 
length L is implicitly computed from 


L ~ V e /x for 2 - D flow 


( 15 ) 
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whore 



(16) 


The non-dimensional heating level Q w , defined in 
Ref. 3, includes effects of nonequilibrium boundary- 
layer chemistry and surface catalytic efficiency. 


Test Cases 

Flat Plate 

Test conditions for the flat plate are presented in 
Table 1. The plate is 30 cm long with a sharp leading 
edge radius of 0.1 mm. The flow is laminar, steady, 
and two-dimensional at a = 0. The Mach 6 test case 
is chosen to be compatible with operating conditions 
in the LaRC 20 inch Mach 6 facility for possible fu- 
ture testing. The Mach 4 and Mach 2 test conditions 
with ReL =2.7 10 6 match the total enthalpy of the 
Mach 6 test condition. The sensitivity of the results to 
Reynolds number was checked by introducing a varia- 
tion in freestream density at Mach 2 and Mach 4. The 
lower Reynolds number also provides a better match 
to local test conditions for the RLV test case discussed 
in the next subsection. The high Reynolds number, 
Mach 2 case most closely matches the conditions as- 
sumed in the triple-deck analysis. 

The wall temperature T~ is set to 300 K. Three 
jumps in surface temperature acros xj, corresponding 
to A T w /T~ equal to -0.1, 0.1, and 0.2 are considered. 
Two values for xj equal to 5.06 cm and 20.455 cm were 
simulated in the Mach 6 case to quantify the influence 
of leading edge effects and boundary-layer growth; the 
pressure and heating relaxation phenomena being in- 
vestigated here associated with the temperature jump 
occur on a comparable scale to these background in- 
fluences. Only the results for the jump at xj = 20.455 
cm are presented here, and values of Rei presented in 
Table 1 are based on this location. Comparisons of nu- 
merical simulation and triple deck theory at x j = 5.06 
are qualitatively and quantitatively similar to the com- 
parisons that follow at xj = 20.455 cm. 


Table 1 - Flat Plate Test Conditions 


Moo 

ReL 

Poo 

Uoo 

Too 

r T 

2 

2.7 10 B 

3.433 KT 1 

654.6 

265.9 

300. 

2 

2.7 10 4 5 

3.433 10" 2 

654.6 

265.9 

300. 

4 

2.7 10 6 

1.263 10 -1 

881.6 

120.4 

300. 

4 

2.7 10 5 

1.263 10 -2 

881.6 

120.4 

300. 

6 

2.7 10 6 

5.850 10 -2 

956.0 

63.0 

300. 


Reusable Launch Vehicle (RLV) Test Case 

The Reusable Launch Vehicle (RLV) is the B1001 
Lockheed configuration defined in Phase I design. 
Only the windside forebody is considered in this test. 
The freestream test conditions, originally considered 
in the catalytic jump paper, 3 are defined in Table 2. 
The flow simulation is at a = 45 deg. and 79.6 km 


altitude. Surface temperature is held fixed at 900 K in 
the baseline CFD simulation. A five-species model for 
air and a fully catalytic wall boundary condition are 
applied. 

A temperature jump of 90 K and 900 K ( &T W {T W = 
0.1 and 1.0, repectively) is imposed at an axial loca- 
tion 25.3 m downstream of the nose. This location 
occurs near the peak of the recompression following 
the overexpansion of the flow from the nose along the 
centerline at a = 45 deg. The imposed wall tempera- 
ture of 900 K is the approximate value of the radiative 
equilibrium wall temperature in the vicinity of the 
jump. The parameter A T w /T w — 1.0 is not small 
in the context of triple-deck approximations; however, 
(A T w /T w )(T w /Tad) is small in this application and 
so the theory is tested in this range. 

Two locations across the temperature jump, the first 
along the windside centerline and the second at an 
outboard location just inside of the chine, are stud- 
ied. Boundary-layer edge conditions at each location 
are recorded in Table 2. The equivalent length derived 
from the local baseline heating rate at these locations 
and integral boundary-layer theory as discussed previ- 
ously is also recorded in Table 2. It is this equivalent 
length that is used to define Rej> 

Table 2 - RLV Local Test Conditions 



Numerical Simulation 

Code Description 

The computational fluid dynamic (CFD) solutions 
are provided by the Langley Aerothermodynamic 
Upwind Relaxation Algorithm (LAURA). 5-7 Com- 
parisons to experimental data for hypersonic flows 
in air are documented in the literature. 8-12 The 
code employs upwind-biased, point- and line-implicit 
relaxation. Inviscid fluxes are approximated with 
Roe’s averaging, 13 eigenvalue limiting (similar to 
Harten 14 ), and Yee’s symmetric total variation dimin- 
ishing scheme. 15 Viscous fluxes are approximated with 
central differences. 

Although a space marching code would normally be 
sufficient for describing flow on a flat plate it is not 
appropriate for the present study because identifica- 
tion of upstream influences of the surface temperature 
jump is of primary concern. Consequently, global re- 
laxation of the governing equations is employed. 

Grid Refinement - Flat Plate 

Sequences of successively finer streamwise grids were 
used to resolve the flow in the vicinity of xj. The ini- 
tial grid used 200 cells from the leading to trailing 
edge. Subsequent grids picked up converged solutions 
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0.08 H 



x, m 


Fig. 3 Test case grid for simulation with T w = 300. 
K and factor 100 increase in streamwise resolution 
around xj = 20.455 cm in embedded block. 

of the flow sufficiently upstream of xj to preclude 
any upstream effects of the jump on the inflow pro- 
files. The simulations used 128 cells across the shock 
layer; comparisons to a 64 cell simulation indicates 
grid convergence. The outer boundary is aligned with 
the shock emanating from the leading edge. The cell 
Reynolds number at the wall averaged 0.25. The maxi- 
mum stretching factor in the normal direction was 1.12 
occuring near the middle of the boundary layer at k = 
50. 

Streamwise resolution on the coarsest grid near xj is 
approximately 3 mm. The next grid in sequence uses 
a single cell deep feeder block and a single, refined 
block (200 x 128) extracted from the coarser block 
to resolve upstream and downstream influence of the 
temperature jump at xj. It has streamwise resolution 
near xj of approximately 0.3 mm. 

The finest grid uses a single cell deep feeder block 
and 3 active blocks (89 x 128), (128 x 128), and (90 
x 128) resolving the shock layer near xj. The sec- 
ond of these active blocks provides another factor 10 
increase in streamwise resolution near xj to approxi- 
mately 0.03mm. An overview of the fine grid is shown 
in Fig. 3. 

Grid Refinement - RLV 

Streamwise grid refinement progressed in the same 
manner as the flat plate. The final computational do- 
main includes a windside section of the RLV extending 
from 10 m to 23 m behind the nose. The final surface 
grid and a representative surface heating distribution 
are presented in Fig. 4. The inflow boundary is out- 
side the influence of the temperature jump. The fine, 
streamwise mesh spacing across the jump is approx- 
imately .008 m so that the local Reynolds number 
based on streamwise mesh spacing is approximately 



x, m 


Fig. 4 Heating and surface mesh on RLV test case 
with A T w /T w = 1.0. 

70. A corresponding value for Reynolds number based 
on streamwise mesh spacing on the flat plate is ap- 
proximately 40 for the low Reynolds number tests. 

The grid refinement process referenced an available 
surface grid of the RLV, not the original CAD de- 
scription. Consequently, the grid refinement produces 
a faceted surface in regions where the curvature is 
nonzero in the streamwise direction. Faceting is most 
evident near the outboard chines when viewing surface 
pressure contours. Nevertheless, the nondimensional- 
ization used herein subtracts out the faceted baseline 
reference results for pressure and heating. 

Results and Discussion 

Pressure Field above the Surface 

Pressure contours across the boundary layer of the 
flat plate at = 6 and ReL — 2.7 10 5 6 are presented 
in Fig. 5. The constant wall temperature case is in 
Fig. 5a and the temperature jump at xj = 20.455 
cm is in Fig. 5b. Vertical lines in this figure indicate 
computational block boundaries of the fine grid region 
surrounding the discontinuity. 

The background pressure field in Fig. 5a is nearly 
constant and a very weak oscillation across the bound- 
ary layer edge manifests itself as a jagged profile. If we 
take a cut across the boundary layer and weak shock 
emanating from the leading edge there is no indica- 
tion of a pressure oscillation in Fig. 6a. However, in 
the highly expanded scale of Fig. 6b an oscillation is 
visible at the boundary layer edge. Numerical exper- 
iments on grid stretching and limiters failed to reveal 
a numerical cause for this profile; we note only that 
relative mimima and maxima follow streamlines along 
the boundary layer edge to the highly curved portion 
of the leading edge shock. Whatever their cause, the 
magnitude of the oscillations is two orders smaller than 
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a) without temperature jump 



b) with temperature jump 

Fig. 5 Comparison of pressure contours around 
xj = 20.455 cm in embedded block with and without 
temperature jump. 

the magnitude of pressure perturbations caused by the 
temperature jump. 

The imposition of an abrupt 10% increase in wall 
temperature at xj induces changes in pressure and 
heating that are felt upstream and in the wake of the 
discontinuity. In Fig. 5b, the pressure jump is largest 
at xj and the perturbation moves up to the boundary- 
layer edge. Coalescing pressure waves followed by an 
expansion are observed trailing off to the right. The 
pressure oscillation from the simple, constant wall tem- 
perature boundary condition disappears, masked by 
the larger gradients in induced pressure. 

Pressure contours in the symmetry plane across the 
RLV shock layer associated with temperature jumps 

6 



a) profile across shock layer 



b) profile across boundary layer 

Fig. 6 Pressure profiles across layer highlighting 
magnitude of oscillation near the boundary-layer 
edge. 

A T w /T w equal to 0.0, 0.1, and 1.0 are presented in Fig. 
7. The strengthening compression and expansion with 
increasing A T w are evident in this figure. A reflected 
wave off the bow shock strikes the surface at x — 20 
m. for the highest temperature jump. The upstream 
influence of the discontinuity in all cases disappears 
within .5 m.. 

Comparison of Triple-Deck Theory and CFD 

Flat Plate 

A direct comparison of the present triple-deck the- 
ory with the Navier-Stokes numerical solution for the 
non-dimensional scaled pressure distribution p(x) is 
shown in Fig. 8 for the M& = 2, high Reynolds num- 
ber case. A range of positive and negative temperature 

ll 
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Fig. 7 Pressure contours in symmetry plane asso- 
ciated with different temperature jumps at xj for 
R.LV test case. 



Fig. 8 Effect of temperature jump amplitude on 
interactive pressure distribution for Moo = 2 and 
Re L = 2,700,000. 

jumps are considered. It is seen that the shape, ex- 
tent and magnitude of the pressure perturbation are 
all very accurately predicted by the triple-deck theory 
throughout both the upstream influence and down- 
stream wake regions of the interaction. The same 
excellent comparisons are also obtained for the cor- 
responding non-dimensional heat transfer predictions 
q as shown in Fig. 9. The theory is clearly able to 
properly capture the important sharp peaking of the 
local heating that occurs immediately downstream of 
the temperature jump. The slow subsequent approach 
in the wake toward the proper asymptotic value per- 
taining to the new wall temperature is also captured. 

The application of the leading asymptotic approxi- 

7 


M. = 2 Re_ xl = 2,700,000 



M_ = 2 Re_ x| = 2,700,000 



Fig. 9 Effect of temperature jump amplitude on 
interactive heating distribution for M x = 2 and 
ReL = 2, 700, 000. 

mation (e — > 0) of triple-deck theory at finite Reynolds 
numbers introduces some error that may be assessed 
by comparing with the Navier-Stokes numerical solu- 
tions. This comparison is presented in Fig. 10, where 
the p(x ) distributions for the = 2 test case with 
A T w /T w = 0.1 at two different Reynold numbers is 
compared to triple-deck theory. Note that all Reynolds 
numbers collapse to a single curve in the triple-deck 
theory, which is most accurate for ReL -t oo. It can 
be seen that the finite Reynolds number effect serves 
to slightly reduce the streamwise extent of the inter- 
action in x primarily in the upstream influence region, 
the reduction increasing with decreasing ReL. This 
effect stems from the fact that in contrast to what is 
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Triple Deck (Re L -> «■) 



Fig. 10 Finite Reynolds number effect on interac- 
tive pressure field. 

assumed in the triple-deck model, the upstream re- 
gion of any real interaction depends on conditions at 
the upstream start of the interaction rather than at 
X=L; as discussed by Stewartson, this leads to a mod- 
est Reynolds number effect that reduces the effective 
upstream influence distance in x. The present results 
for the temperature jump problem are also in qual- 
itative agreement with an earlier theoretical study 16 
of finite Reynolds number effects on oblique shock- 
generated interactions. 

Turning to a consideration of the Mach number ef- 
fect, Fig. 11 presents a comparison of the theoretical 
p(x) distribution (which is independent of M^) with 
the CFD solutions obtained at Mach numbers of 2, 4, 
and 6 . It is seen that as the freestream Mach num- 
ber increases and enters the hypersonic regime, the 
numerical simulation of the exact equations predicts 
a progressively more significant reduction in upstream 
influence along with a modest reduction in the corre- 
sponding maximum pressure. Analysis indicates that 
this Mach number effect is in fact the combined result 
of nonlinear effects in both the inviscid pressure-flow 
deflection relationship and the temperature jump in- 
teraction (neither of which are accounted for in the 
present linearized treatment) . Interestingly there is lit- 
tle attendant Mach number effect on the downstream 
wake behavior. 

The corresponding effects of Mach number on heat- 
ing across the temperature jump are shown in Fig. 12. 
As Mach number increases the upstream influence of 
the temperature jump on heating decreases in x coor- 
dinates faster than predicted by theory. This trend is 
the same as observed in Fig. 11 for the upstream influ- 
ence of Mach number to p(x). A magnified view of the 
downstream asymptotic behavior of q(x) as a function 
of Mach number is shown in Fig. 12b. The down- 



Fig. 11 Mach number effect on interactive pres- 
sure field. 

stream behavior of q(x ) is not independent of in 
the triple-deck theory because of the functional depen- 
dence of C 5 on 0 and Cref ■ All three Mach number 
simulations show nearly a 20 % difference in q at the 
most downstream x location. This trend will be evi- 
dent in the RLV simulations as well. The magnitude of 
the parameter C 5 has the most influence on this com- 
parison in this range where ex 1 / 3 is of order 1. Note 
that the limit for q{x) as f — > 00 is correct in Eq. 10. 

RLV 

The RLV test case is intended to show how unmod- 
ified triple-deck theory performs in a regime that is of 
practical interest to the thermal protection system de- 
sign community. The RLV simulation includes effects 
of chemical nonequilibrium (5 species air model) and a 
fully catalytic wall boundary condition. The limiting 
form of Eq. 10 does not account for species diffusion to 
the surface. The results obtained here will be used to 
guide more comprehensive analytic modeling of both 
temperature and catalytic jump problems in the con- 
text of triple-deck theory. In spite of these caveats, the 
theory captures much of the interaction quite well. 

The strong three-dimensionality of the flow is evi- 
dent in the surface heating contours of Figs. 4 and 13. 
The extent of the interaction can be discerned by com- 
paring the heating rates for A T w /T w equal to 0 and 
0.1 in Fig. 13. It is also clear that this perturbation in 
temperature is very low relative to the adiabatic wall 
temperature for this case - thus evolving an extended 
study with a larger perturbation A T w /T w = 1.0. 

A comparison of triple-deck theory and numerical 
simulation for the pressure and heating across the tem- 
perature jump A T w /T w = 0.1 is presented in Fig. 
14 for the centerline location. The agreement in up- 
stream influence and downstream asymptotic behavior 
of pressure is very good. Only the peak level in pres- 
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Fig. 12 Mach number effect on interactive heating 
distribution for AT W /T W =0.1 and ReL = 2,700,000 
on expanded x scale. 


sure disturbance at the jump is overpredicted, possibly 
because the theory does not account for the varying 
heat capacity across the boundary layer. The up- 
stream heating influence is quite similar to the agree- 
ment found at Mach 2 for the flat plate. The relaxation 
to the asymptotic value of heating downstream of the 
jump appears offset - though in magnitude there is less 
than a 2% decrement in the asymptotic heating rate 
associated with the 10% increase in surface tempera- 
ture. 

A comparison of triple-deck theory and numerical 
simulation for the pressure and heating across the tem- 
perature jump A T w /T w = 0.1 is presented in Fig. 
15 for the outboard location. The flow in this re- 



Fig. 13 Surface heating rates for RLV test case 
with and without temperature jump, AT W /T W = 0.1. 

gion is more strongly three-dimensional and the local 
Reynolds number is a factor 2.5 smaller than at the 
centerline. The comparison follows the same trends as 
at the centerline but the magnitude of differences is 
larger. 

A comparison of triple-deck theory and numerical 
simulation for the pressure and heating across the 
larger temperature jump A T w /T w = 1.0 is presented 
in Fig. 16 for the centerline location. As noted pre- 
viously, this temperature perturbation can no longer 
be considered small in the context of the triple-deck 
formulation. Nevertheless, it is small relative to the 
adiabatic wall temperature. The pressure distribution 
scales exactly with A T w /T w - the corresponding curves 
for theory and CFD in Fig. 14 are almost identical 
to those in Fig. 16. The CFD heating distribution 
is approximately 10% higher than theory at x = 20, 
a location near the downstream edge of the fine grid 
in Fig. 4. The CFD heating distribution is approxi- 
mately 6% higher than theory at x — 200, a location 
near the outflow boundary of the larger domain. 

In all cases, including the flat plate, the CFD sim- 
ulation approaches its asymptotic level faster than 
predicted by triple deck theory in the range 1 < x < 
Order(lOO). The behavior of q in this domain is domi- 
nated by the term in Eq. 10 which is of Order(l) 
in this range. 

Concluding Remarks 

This paper has presented a new analytical theory of 
wall temperature jump-generated interactions in su- 
personic flow that captures all the major physics of 
the resulting pressure and local heat transfer behav- 
ior, including upstream influence. The long range goal 
of this work is to provide analytic tools coupled with 
computational data sets to be used in the design of 
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Fig. 14 Triple-deck theory applied to RLV test 
case with AT W /T W = 0.1 on the windside plane of 
symmetry. 

hypersonic vehicles in an Intelligent Synthesis Envi- 
ronment (ISE). Given a single, baseline CFD solution, 
the analytic tools will allow a designer to investigate 
effects of continuous global changes or discontinuous 
local changes to thermal protection system surface 
properties on heating. 

Comparisons of triple-deck theory with numerical 
solutions of the exact Navier-Stokes equations for a 
perfect gas on a flat plate are presented. Good agree- 
ment establishes confidence that the closed form ana- 
lytical results are properly formulated and coded. In 
the case of laminar, high Reynolds number, supersonic 
flow (Moo = 2 and Re^ — 2.7 10 6 ) theoretical pre- 
dictions for the upstream and downstream behavior 

10 




b) heating 

Fig. 15 Triple-deck theory applied to RLV test 
case with AT W /T W = 01 near the outboard location. 

of surface pressure and heating in the vicinity of the 
discontinuity are in excellent agreement with computa- 
tional simulation. Some decrement of theory is noted 
for larger edge Mach numbers (hypersonic) and for the 
asymptotic behavior of heating far downstream from 
the jump. 

A Reusable Launch Vehicle (RLV) test case with an 
imposed surface temperature discontinuity is intended 
to show how unmodified triple-deck theory performs 
in a regime that is of practical interest to the thermal 
protection system design community. The RLV sim- 
ulation includes effects of chemical nonequilibrium (5 
species air model) and a fully catalytic wall boundary 
condition. Although the nonequilibrium physics of this 
problem is not included in the present triple-deck for- 

ll 
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Fig. 16 Triple-deck theory applied to RLV test 
case with AT W /T W = 1.0 on the windside plane of 
symmetry. 

mulation, the theory captures much of the interaction 
quite well. Guidance for adapting a strongly three- 
dimensional, CFD flow simulation (on a grid that is 
too coarse to resolve the interactions across the dis- 
continuity) into the context of a triple-deck analysis is 
provided. 

A problem related to the temperature jump is that 
of dissociated gas flow past a wall with a jump in 
the surface catalycity (heterogeneous recombination 
rate). This situation arises in connection with ther- 
mal protection system (TPS) design for hypervelocity 
vehicles operating in the high-altitude nonequilibrium 
flow chemistry regime where there is a discontinuity in 
the surface material at some station along the body. 


Application of triple-deck theory to this more complex 
problem will be guided by the results obtained herein 
for the RLV application. 
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